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Abstract 

We report a set of measurements of inclusive invariant pr differential cross sections of A*^, a'^, 
S^, and hyperons reconstructed in the central region with pseudorapidity |r7| < 1 and pr up 
to 10 GeV/c. Events are collected with a minimum-bias trigger in pp collisions at a center-of-mass 
energy of 1.96 TeV using the CDF II detector at the Tevatron Collider. As pt increases, the slopes 
of the differential cross sections are similar not only to each other but also to those of mesons, 
which could indicate a universality of the particle production in px The invariant differential cross 
sections are also presented for different charged-particle multiplicity intervals. 
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es containing strange 
ippQ,^ andppU). 



Ever since their discovery in cosmic ray interactions partic 
quarks have been extensively studied at particle colliders (e"'"e~ [3], ep 
The process by which hadrons in general are produced from interactions is an unsolved 
problem in the standard model, and a detailed analysis of production properties of particles 
with different quark flavors and numbers of quarks could pave the way to understanding the 
process from first principles. The data on strange particle production can also be used to 
refine phenomenological models and set parameters, such as the strange quark suppression 
constant in event generators, which have become an integral part of any data analysis. 
Interest in particles containing strange quarks increased with the introduction of the quark- 



d manifest itself as 
tI. To isolate QGP 



gluon plasma (QGP). Formation of quark-gluon plasma in a collision cou 
an enhanced production of strange particles such as kaons and hyperons 
signatures in heavy-ion collision data, understanding the particle production properties from 
simple nucleon interactions is necessary. 

There are ample data on the production ofparticles with one strange quark, but very little 



available on particles wit 



1 two or more 



such as RHIC 



10|, SppS 



m oipc 

m 



Previous studies of hyperons from colliders 



11| , and the Tevatron [12|, |13| were limited by low sample statistics 



and the limited accessible range of hyperon momentum component transverse to the beam 
direction (pt)- In this Letter, we report on a study of the hyperons A° (quark content uds), 
E~ {dss), and Q~ {sss) and their corresponding antiparticles (A°, S"*", and For these 
hyperons, the inclusive invariant pt differential cross sections are measured up to px of 10 
GeV/c, based on ~ 100 million minimum-bias events collected with the CDF II detector. 
The measurements reported here are the current best from any hadron collider experiment 
in terms of statistics and pt range. 



The CDF II detector is described in detail elsewhere [ij]. The components most relevant 
to this analysis are those that comprise the tracking system, which is within a uniform axial 
magnetic field of 1.4T. The inner tracking volume is composed of a system of eight layers of 



silicon microstrip detectors ranging in radius from 1.5 to 28.0 cm 15| in the pseudorapidity 



region jr^l < 2 [l6]. The remainder of the tracking volume is occupied by the Central Outer 
Tracker (COT). The COT is a cylindrical drift chamber containirig 96 sense- wire layers 
grouped in eight alternating superlayers of axial and stereo wires 17|. Its active volume 
covers 40 to 140 cm in radius and \z\ < 155 cm. The transverse- momentum resolution of 
tracks reconstructed using COT hits is o'{pt)/Pt ~ 0.0017/(GeV/c). 
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Events for this analysis are collected with a "minimum-bias" (MB) trigger, which selects 
beam crossings with at least one pp interaction by requiring a timing coincidence for signals in 
both forward and backward gas Cherenkov counters ISj covering the regions 3.7 < |?7| < 4.7. 
The MB trigger is rate-limited to keep the final trigger output at 1 Hz. Primary event vertices 
are identified by the convergence of reconstructed tracks along the beam axis. Events are 
accepted that contain a reconstructed vertex in the fiducial region {z^txl < 60 cm centered 
around the nominal CDF origin {z = 0). When an event has more than one vertex, the 
highest quality vertex, usually the one with the most associated tracks, is selected and it 
is required that there be no other vertices within ±5 cm of this vertex. This selection 
introduces a bias toward high multiplicity events as the instantaneous luminosity increases. 
To combine events collected at different average instantaneous luminosities, we determine a 
per-event weight as a function of the charged-track multiplicity Nch in order to match the 
multiplicity distribution of a data sample where the average number of interactions is less 
than 0.3 per bunch crossing. For the Nch calculation, tracks are required to have a high 
track-fit quality with per degree-of-freedom {x^/dof) less than 2.5, and more than five 
hits in at least two axial and two stereo COT segments. It is further required that tracks 
satisfy \ri\ < 1, impact parameter do [l9| less than 0.25 cm, the distance along the 2;-axis 
(5Zo)between the event vertex and the track position at the point of closest approach to the 
vertex in the r — cj) plane be less than 2 cm, and px > 0.3 GeV/c. The pt selection is to 
minimize the inefficiency of the track-finding algorithm for low momentum tracks. 

We search for A° — )• pn^ decays using tracks with opposite-sign charge and pr > 0.325 
GeV/c that satisfy the x^/dof and COT segment requirements. In this Letter, any reference 
to a specific hyperon state implies the antiparticle state as well. For each two-track combi- 
nation we calculate their intersection coordinate in the r — cf) plane. Once this intersection 
point, referred to as the secondary vertex, is found, the z-coordinate of each track {Zi and 
Z2) is calculated at that point. If the distance |Zi — Z2I is less than 1.5 cm, the tracks are 
considered to originate from a A° candidate decay. The pair is traced back to the vertex and 
we require 5Zq be less than 2 cm, and the do be less than 0.25 cm. To reduce backgrounds 
further, we require the A'' decay length Lao, the distance in the r — cf) plane between the 
primary and secondary vertices, to be greater than 2.5 cm and less than 50 cm. 

The invariant mass Mp^ of the two-track system is calculated by attributing the proton 
mass to the track with the higher p-p, as preferentially expected by the kinematics of a A° 



10 



decay. Figure [T] shows the invariant mass for A° candidates with \ri\ < 1. This distribution 



2 is divided into 23 pt intervals 



20| and the number of A° in each px interval is determined by 



3 fitting the invariant mass distributions using a Gaussian function with three parameters for 

4 the signal and a third-order polynomial for the underlying combinatorial background. The 

5 data in the mass range 1.10 — 1.16 GeV/c^ are fitted. The polynomial fit to the background 

6 is subtracted bin- by-bin from the data entries in the A° mass window (1.111 — 1.121 GeV/c^) 

7 to obtain the number of A" hyperons. This number is divided by the acceptance to obtain 

8 the invariant pt differential cross section as described later. 

9 The fitting procedure is one source of systematic uncertainty. This uncertainty is esti- 

10 mated by separately varying the mass range of the fit, the functional form for the signal 

11 to a double Gaussian, and the background modeling function to a second-order polynomial. 

12 The number of A*^ is recalculated in all pt intervals for each variation. The systematic un- 

13 certainty is determined as the sum in quadrature of the fractional change in the number of 

14 A° from each modified fit. It decreases from ±10% at the lowest pr (1.2 GeV/c) to less 

15 than ±5% for pr > 1.75 GeV/c. 

16 The cascade reconstruction decay mode is — t- A'^tt^ — )■ {pTi^)^^. The previously 

17 reconstructed A° candidates are used, but without the do and 6Zo requirements. We select 
IB A° candidates in the A'' mass window and calculate the coordinate of the intersection point 

19 in the r — (p plane between the A° candidate and a third track. The 2;-axis coordinates at 

20 this point are calculated for the third track (Z3) and the A'' candidate {Z4). The three-track 

21 system is considered a S~ candidate decay if the distance {Z^ — Z^\ < 1.5 cm. We also 

22 require the decay length L^- > 1 cm, and (Lao — L=- ) > 1 cm. Finally, it is required that 

23 the do and 6Zq of the S~ candidate be less than 0.25 cm and 2 cm respectively. 

24 The invariant mass M^o^ is calculated by fixing the mass of the A° candidate to 1.1157 

25 GeV / |2l| and assigning the pion mass to the third track. Figure [1] shows the invariant 

26 mass for H~ candidates with |?7| < 1. 

27 The H~ candidates are divided into 17 pr intervals and the number of S~ in each interval 
2B is determined by fitting the corresponding M\o^ invariant mass distribution using a Gaus- 

29 sian function for the signal and a third-order polynomial for the background. The fitted 

30 background is then subtracted bin-by-bin from the data entries in the signal region (1.31 to 

31 1.33 GeV/c^) to obtain the S~ yield in every pt interval. The systematic uncertainty of the 

32 fit procedure is estimated the same way as for the A° and is found to change by no more 
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than ±5% in all pt intervals. 

To reconstruct Q~ decays we follow the same procedure as for the H~ and apply the same 
selection criteria except that the third track is assigned the kaon mass. The search decay 
mode is — t- J\PK~ — )■ {p7T^)K^. Because of the larger background, the procedure to 
extract the signal yield is slightly different from that in the previous cases. Track pairs 
with Mpjr- in the mass ranges 1.095 — 1.105 and 1.127 — 1.137 GeV/c^ are combined with 
the third track to obtain the invariant mass distribution of the combinatorial background. 
This distribution is subtracted from the M\ok~ distribution after normalizing to the number 
of events in the mass window 1.69 < Mj^o^- < 1.74 GeV/c^. The background subtracted 
M\ok- invariant mass distribution is shown in Fig. [T] 

The distribution is divided into 10 pt intervals, and we use the method described above 
to extract the Q~ signal from the corresponding invariant mass distributions in each pt 
interval within the mass window 1.66 to 1.68 GeV/c^. The systematic uncertainty due to 
the fitting procedure is also calculated in a similar manner as S~, with the exception of 
using a double Gaussian variation because of low statistics. The overall uncertainties 
are about ±10% for all Pt intervals. 

The geometric and kinematic acceptance is estimated with Monte Carlo (MC) simula- 
tions. The MC data of a resonance state are generated with fixed pt corresponding to 14 



points 



201 ] ranging from 0.75 to 10 GeV/c and flat in rapidity \y\ < 2. A generated resonance 



is combined with either one or four non-diffractive inelastic MB events generated with the 
PYTHIA ^1 generator. Although the average number of interactions in our data sample is 
a little less than two, the default acceptance is calculated from the MC sample with four 
MB events and the difference of the acceptance values between the two samples is one of 
our systematic uncertainties. Based on a study with tracks from Kg decay, the sample with 
four MB events reproduces the low px tracking efficiency in data well within the systematic 
uncertainty 

The detector response to particles produced in the simulation is modeled with the CDF II 
detector simulation that in turn is based on the geant-3 MC program 23[. Simulated events 
are processed and selected with the same analysis code used for the data. The acceptance 
is defined as the ratio of the number of reconstructed resonances with the input pt over 
the generated number, including the branching ratio. Acceptance values are calculated 
separately for the particles and their corresponding antiparticles and the average of the 
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1 two is used as the default value, since the acceptances for the two states are similar. The 

2 acceptance values obtained for the 14 pr points are fitted with a fourth order polynomial 

3 function and the fitted curve is used to correct the numbers of each hyperon state in the 

4 data. 

5 The modeling of the MB events overlapping with the examined resonance and the selec- 

6 tion criteria applied contribute as a systematic uncertainty to the acceptance calculation. 

7 The contribution from the former has already been mentioned. Acceptance uncertainties 

8 due to the selection criteria are studied by changing the selection values of the variables used 

9 to reconstruct the resonances. The variables examined are pr, |Zi — Z2I, IZ3 — Z4I, SZq, do 

10 and the decay lengths. For each variable other than Pt, two values around the default value 

11 are typically chosen. One value is such that it has little effect on the signal, and the other 

12 reduces the signal by ~20 to 30%. The default pr selection value is 0.325 GeV/c, and it is 

13 changed to 0.3 GeV/c and to 0.35 GeV/c. 

14 For each considered variation, a new acceptance curve and the number of resonances as a 

15 function of pr are obtained, and the percentage change between the new pr distribution and 

16 the one with the default selection requirements is taken as the uncertainty in the acceptance 

17 for the specific pr interval. The square root of the quadratic sum of the uncertainties from 

18 each variation is taken as the total conservative uncertainty on the acceptance in a given 

19 Pt bin. The systematic uncertainty associated with the Q~ hyperon acceptance is derived 

20 from the E~ uncertainty estimate since the reconstruction follows the same criteria. This 

21 acceptance uncertainty is added quadratically to the systematic uncertainty due to the fitting 

22 procedure, described earlier, to give the total systematic uncertainty. 

23 For the A*' case, the acceptance uncertainty decreases from about 25% at ~ 1 GeV/c 

24 to 10% at Pt ~ 2 GeV/c and then rises again slowly to 15% for pt > 7 GeV/c. The 

25 corresponding acceptance uncertainty for the H~ {^~) case decreases from about 15% (20%) 

26 at Pt ~ 2 GeV/c to 10% (15%) for pt > 4 GeV/c. 

27 The inclusive invariant pr differential cross section for each hyperon resonance is calcu- 

28 lated as Ed^a/dp'^ = {amb/ N event) d^ N / ApTdpTdydcj) = {amb/'27iNeyent)^N/ApTApTAy where 

29 (Tmf, is our MB trigger cross section, Nevent is the number of weighted events, AA^ is the num- 

30 ber of hyperons observed in each pt interval { Apj-) after background subtraction, A is the 

31 acceptance in the specific pr interval, and Ay is the rapidity range used in the acceptance 

32 calculation (-2 to 2). 



13 



TABLE I: The results of power law function fits to the px differential cross sections described in 
the text and shown in Fig. [2] for -pr > 2 GeV/c. The parameter pQ is fixed to 1.3 GeV/c in all fits. 
The values are from Ref. 25| at ^/s = 1.8 TeV. The uncertainties shown do not include the 
MB cross section uncertainty [24] . The last line of the table gives the per degree-of- freedom of 



the fit to data. 



Parameter (units) K^g [25] A° 



^ (mb/GeV^c^) 45 ± 9 210 ± 25 14.9 ± 2.5 1.50 ± 0.75 
Po (GeV/c) 1.3 1.3 1.3 1.3 

n 7.7 ± 0.2 8.81 ± 0.08 8.26 ± 0.12 8.06 ± 0.34 
X^/dof 8.1/11 5.7/15 15.8/15 10.5/7 



Figure |2] shows the results for the pt differential cross section for the three hyperon 
resonances. The uncertainties shown for each data point include the statistical and all 



systematic uncertainties described above, except the one associated with amb [24 1. 

The Pt differential cross section is modeled by a power law function, Aiv n)'^ J (pr + Po)"", 
for Pt > 2 GeV/c. In order to compare with the previous CDF Kg result 5|, l25|], po is fixed 



at 1.3 GeV/c, and the results are shown in Tab. [B The data below pj- ~ 2 GeV/c cannot 
be described well by the power law function even if pq is allowed to float. For this region, 
the data are better described by an exponential function, B exp[—b ■ Pt]- The results of 
this fit are shown in Tab. [Tll and the slope 6 of A is consistent with previous measurements 

nn ^ p 

[121 . 1131 ]. The b values depend on the range of the fit but are about two, which corresponds 
to an average pt of 1 GeV/c under the assumption that the fit can be extrapolated down to 
p^ = GeV/c. 

The plots on the right side of Fig. [2] show the ratio of the pt differential cross sections 
for and A", and Q~ and A". In the H^/A*^ ratio there is a rise at low pt, and the ratio 
reaches a plateau at pr > 4 GeV/c. It should be noted that the A° cross section also includes 
A° production from the decay of other hyperon states — )■ A°7, S"^, S° and ^). Due to 
the short S° lifetime, A° from S° decays cannot be separated from direct A° production. 
Simulations of cascade decays indicate that ~ 50% of A° from S decays will satisfy our A° 
selection criteria, with the fraction of A" fairly independent of S pt- The ratio plots in Fig. [2] 
are fitted to a constant, and the value 0.17 ± 0.01 is obtained for S^/A° and 0.025 ± 0.002 
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TABLE II: The results of exponential function fits to the px differential cross sections shown in 
Fig. [2] for the pT ranges given in the second row. The uncertainties shown do not include the MB 
cross section uncertainty [24]. The last line of the table gives the P^r degree-of-freedom of the 
fit to data. 

Parameter (units) A° A° Q.^ 

PT range (GeV/c) [1.2, 2.5] [1.2, 4] [1.5, 4] [2, 4] 

B (mb/GeV^c^) 4.68 ± 1.04 3.16 ± 0.35 0.16 ± 0.04 0.024 ± 0.011 
b (GeV'^c) 2. 30 ±0.12 2. 10 ±0.04 1.75 ±0.08 1.80 ±0.19 
X^/dof 1.0/7 7.2/12 4.0/8 6.3/3 



for n-/A^. 

The plots in Fig. |2] clearly show that the cross sections depend on the number of strange 
quarks. However, the similarity of the n values in Tab. [U along with the measured value 
n = 8.28 ± 0.02 for all charged particles [26], indicate that the slope of the pt differential 
cross sections are similar in the high pt region. This could be an indication of a universality 
in particle production as Pt increases [23]. This is in contrast to the low pt region where 
the slope exhibits a strong particle type dependence 



Figure [3] shows the pt differential cross sections for two charged-particle multiplicity 
regions, Nch < 10 and Nch > 24. N^h = 24 (10) corresponds to dN/di] ~ 16 (7), corrected 
for the track reconstruction efficiency and unreconstructed tracks with pt < 0.3 GeV/c |26 |. 
Due to the low sample statistics, distributions are only shown for A° and We 
observe a correlation between high pt particles and high multiplicity events. This is a 
general characteristic independent of the particle types. 

In summary, the production properties of A°, H~, and Q~ hyperons reconstructed from 
minimum-bias events at -y/i = 1.96 TeV are studied. The inclusive invariant px differential 
cross sections are well modeled by a power law function above 2 GeV/c pr- With fixed po, 
the fit parameter n decreases from 8.81 ± 0.08 (A°) to 8.06 ± 0.34 (f2^). The low pt regions 
are modeled by an exponential function. The exponential slope, b, decreases by ~ 15% from 
A° to Q^. The cross section ratios S~/A° and Q~/A^ are presented as a function of pt- 
Although the ratios exhibit a strong dependence on the number of strange quarks, the n 
values of the hyperons. Kg and all charged particles are within ~ 10% of each other. This 
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could be an indication that the production process which determines the pt of these particles 
depends little on the particle type as pt increases. We also find the hyperon px differential 
cross sections fall off faster with px for low multiplicity events than for high multiplicity 
events. 



xlO' 




(GeV/c') M^.^ (GeV/c') M^.^ (GeV/c') 

FIG. 1: Reconstructed invariant mass distributions for Mp^ (left), M^o^ (center), and Mj^oj^ 
(right). The background has been subtracted from the M/^oj^ distribution. The solid hues are 
fitted curves, a third-degree polynomial for the background and either a double (Mp^ and Mj^o^) 
or single {M^ok) Gaussian function to model the peak. The widths reflect the tracking resolution 
and are consistent with the widths from MC simulation. 
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FIG. 2: The pt differential cross sections for A'^, and fl^ within |r/| < 1 (left). The solid curves 
are from fits to a power law function, with the fitted parameters given in Tab. [H The ratios of 
S~/A^ and Q^/A^ as a function oi pT (right). 
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FIG. 3: The pT differential cross sections for two charged-particle multiplicity regions, N^h < 10 
and Nch > 24. Distributions for A'^ are shown on the left while distributions for are shown on 
the right. 
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